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Kis semileptonic form factor from 2 + 1 flavor lattice QCD 
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We present the first results for the Ki^ form factor from simulations with 2+1 flavors of dynamical 
domain wall quarks. Combining our result, namely /+(0) = 0.964(5), with the latest experimental 
results for Ki^ decays leads to \Vus\ = 0.2249(14), reducing the uncertaintity in this important 
parameter. For the 0(p^) term in the chiral expansion we obtain A/ = —0.013(5). 

PACS numbers: 11.15.Ha,12.15.Hh,12.38.Aw,12.38.-t,12.38.Gc,13.20.Eb 



The increasing precision with which the unitarity of 
the CKM matrix [l| can be tested is an important tool 
for exploring the limits of the Standard Model. One such 
unitarity relation is 



IK, 



IK. 



|' + |Kb|' = l 



(1) 



whose uncertainty is dominated by the precision of | Ks I ■ 
In order to obtain |Ks| from experimental measurements 
of the rate for an s — > u decay process, it is necessary 
to quantify the corresponding non-perturbative QCD ef- 
fects. In this paper we present first lattice results from 
simulations with 2+1 flavors of domain wall quarks, 
which respect chiral and flavor symmetries to high ac- 
curacy, for the evaluation of the form factor /+(0) nec- 
essary to determine |Ks| from K t\Ivi (Ktz) semilep- 
tonic decays. Precise knowledge of /+(0) is crucial also 
for deducing |Vtd| from a measurement oi K ^ tt^vv. 

Our determination of /+(0) includes estimates of all 
systematic errors (chiral and extrapolations, dis- 
cretization and finite volume effects), and reduces the 
combined theoretical and experimental error in | Ks | from 
the PDG(2006) result of 0.2257(21) to ^ 



iKsl = 0.2249(14) 



(2) 



The combination |Ks/+(0)| can be obtained from the 
experimental rate for Ki^ decays 



^5(7(2), Sew^ Asm contain the isospin breaking, short 
distance electroweak and long distance electromagnetic 
corrections, respectively. /+(0) is the form factor defined 
from the K ~f n matrix element of the weak vector cur- 
rent, Vf^ — sj^u, evaluated at zero momentum transfer 

{7rip')\V^\Kip)) = {j>,+p'^)f+{q^) + (p^ -p;)/_(g2) , 



where = {p - p'f . PDG(2006) quotes 
|Ks/+(0)| =0.2169(9) , 



(4) 
(5) 



hence in order to obtain |Ks| at a precision commen- 
surate with current experiments, we need to determine 
/+(0) with an error of less than 1%. 

In chiral perturbation theory (ChPT), /+(0) is ex- 
panded in terms of the light pseudoscalar meson masses 

/+(0) = l + /2+/4 + ..., {fn=0{ml^^^)) . (6) 

Current conservation ensures that in the S'J7(3)flavor limit 
— li hence /2 and /4 are small. Additionally, as a 
result of the AdemoUo-Gatto Theorem which states 
that /2 receives no contribution from local operators ap- 
pearing in the effective theory, /2 is determined unam- 
biguously in terms of m,r, mK and /,r, and takes the 
value /2 = —0.023 at the physical values of the meson 
masses @ . Our task is now reduced to one of finding 



K' 



1927r3 



IS 



EW 



(3) 



X [1 + 2Asc,(2) + 2A£m] |K.n/+(0)P , 

where / is the phase space integral which can be evalu- 
ated from the shape of the experimental form factor, and 



A/ = /+(0)-(l + /2) 



(7) 



Until recently, the canonical estimate of A/ = —0.016(8) 
was due to Leutwyler & Roos (LR) [6], whereas more 
recent ChPT based phcnomenological analyses favor a 
value consistent with zero (see Table |l|. These deter- 
minations, however, require model input; the 50% error 



For a recent review of lattice determinations of /x (necessary 
for the determination of \Vus\ from K£2 leptonic decays) and 

previous computations of /+(0) see 0|. ^ A more recent analysis finds |y„s/+(0)| = 0.21673(46) 
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TABLE I: Summary of ChPT and lattice results. + Results 
in conference proceedings only. * Used slope of experimental 
form factor as input. ^ Information not provided. 



Ref. 


/+(0) 


A,/ 


[GeV] 


a [fm] 


Nf 


M 


0.961(8) 


-0.016(8) 








m 


0.978(10) 


+0.001(10) 








[8] 


0.984(12) 


+0.007(12) 








[9] 


0.974(11) 


-0.003(11) 








[10] 


0.960(5)(7) 


-0.017(5)(7) 


> 0.5 


0.07 





[13] 


0.968(9)(6) 


-0.009(9) (6) 


> 0.49 


0.12 


2 


[l^t. 


0.962(6)(9) 


-0.015(6)(9) 


t 


t 


2+1 


[12]t 


0.967(6) 


-0.010(6) 


> 0.55 


0.09 


2 


[l^t 


0.965(2) 


-0.012(2) 


> 0.5 


0.08 


2 


This work 


0.964(5) 


-0.013(5) 


> 0.33 


0.114 2+1 



in the LR result, for example, was estimated within the 
context of a simple quark model. Hence a model indepen- 
dent determination of A/ with a reliable error estimate 
is necessary. We compile recent lattice and phenomeno- 
logical results in TableHl Our lattice calculation has been 
discussed in preliminary form in [15] and we now finalize 
our results with the inclusion of the complete set of data 
and a careful estimate of all systematic errors. 

We simulate with Nf = 2 + 1 dynamical flavors gen- 
erated with the Iwasaki gauge action [l^ at P — 2.13, 
which corresponds to an inverse lattice spacing a^^ = 
1.73(3) GeV (a = 0.114(2)fai) [13, El, and the do- 
main wall fermion action |19j with a residual mass of 
awres = 0.00315(2) The simulated strange quark 

mass, arus = 0.04, is close to its physical value [l8|, and 
we choose four values for the light quark masses, am„d, 
which correspond to pion masses as light as 329 MeV 
[13, Ell. The calculations are performed on two volumes, 
16^ ((1.83)3 fm^) and 24^ ((2.74)^ fm^), at each quark 
mass, except the lightest mass which is only simulated 
on the larger volume. Simulation details are summarized 
in Table [Til and more details can be found in 17, isj . 

We start by rewriting the vector form factors given in 
(HI to define the scalar form factor 



(8) 



which can be obtained on the lattice at qL^^ — {rriK — 
m-„Y ■^ith high statistical accuracy 13, 2^- In Table Hill 
we present our results for /o(9max) fo'" each of the simu- 
lated quark masses and volumes. 

For each quark mass, in addition to evaluating /o('Z^) 
at — g^ax: ^e determine the form factor at several 
negative values of q^, allowing us to interpolate to — 0. 
Specifically, in the notation of (jl]), we evaluate the form 
factor with \p'\ — 0, |p| = pl or \p\ = \plpL where 
PL = 27r/L and L is the spatial extent of the lattice, and 



TABLE II: Simulation parameters: bare light quark mass 
(am„d), pion (mir) and kaon {ttik) masses for both volumes. 





16^ 


X 32 


24^ 


X 64 


aruud 


[GeV] 


niK [GeV] 


[GeV] 


niK [GeV] 


0.03 


0.674(11) 


0.723(12) 


0.671(11) 


0.719(12) 


0.02 


0.557(9) 


0.666(11) 


0.556(9) 


0.663(11) 


0.01 


0.428(7) 


0.614(10) 


0.416(7) 


0.604(10) 


0.005 






0.329(5) 


0.575(9) 



also with \p \ — 0, \p'\ — pl or \p'\ — V^pl ■ To obtain 
the /o('Z^) we use standard ratio techniques 13, 13, 20| . 
which do not require normalization of the vector current. 

In order to gain the maximum amount of information 
from limited data, we perform a simultaneous fit to both 
the q^ and quark mass dependencies using the ansatz 



fo{q 



2 2 



l + f2 + {m%- ml)^{Ao + Ai{ml + ml)) 



1 - qy{Mo 



K 



m2))2 



(9) 



with four fit parameters vIq, ^i, -^o, -^i, a-nd where /2 
is the NLO term appearing in the chiral expansion ([6]), 
evaluated by inserting the lattice results for TOtt , rriK and 
the physical value for (132 MeV) into the expression 
appearing in ChPT [3] at each quark mass 

The expression ([9]) is well motivated since we know 
from the Ademollo-Gatto Theorem that to leading order 
A/ cx (m^ — m„d)^, hence we expect 

/o(0) = l + f2 + {m\-ml)\Ao + A^{m\ + ml)), (10) 

which incorporates the correct >5'C/(3)flavor limit, /+(0) = 
1, to be a good phenomenological ansatz for the mass 



TABLE III: Results for /o(gmax) where g,nax ~ {t^k 





16^ 


X 32 


24^ 


X 64 


airiud 


gmax [GeV^] 


/o('7max) 


*nax [GeV^] 


/o(9max) 


0.03 


0.00233(4) 


1.00035(3) 


0.00235(4) 


1.00029(6) 


0.02 


0.01178(24) 


1.00241(19) 


0.01152(20) 


1.00192(34) 


0.01 


0.03475(66) 


1.01436(81) 


0.03524(62) 


1.00887(89) 


0.005 






0.06070(107) 


1.02143(132) 



Here we note that by using fn in /2 , we are following convention 
Q. The true SU(3) LEG /o, however, is Ukely to be somewhat 
smaller, resulting in a larger and more dominant contribution 
coming from the corresponding /2, and hence a more appar- 
ent convergence of the ChPT. Our lattice calculation determines 
1 — /+ directly and will only differ slightly in the extrapolation 
ansatz. 
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FIG. 1: Scalar form factor, /o(g^), together with the simulta- 
neous fit of © as described in the text. 
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FIG. 2: Scalar form factor, /o(0), together with the simulta- 
neous fit (solid line) on the 24"^ data (red circles) using ([9}. 



dependence of /o(0) = /+(0). This motivates the numer- 
ator in ^ , while the denominator comes from simply in- 
cluding a quark mass dependence into the standard pole 
dominance form 



/o(g') = /o(0)/(l-gVM2) 



(11) 



where M is a pole mass, which has been shown to de- 
scribe the g^-dependence of lattice results of /o(<Z^) very 
well [HEl. 

The traditional approach of sequentially interpolating 
in PT|) followed by chiral extrapolation of /+(0) fTU]) 
should agree with our simultaneous fit Fitting the 
24'^ data only yields excellent agreement (shown in the 
final two rows of Table HV)) . with a reduced error evident 
in the simultaneous fit, which we therefore take as our 
best result. For the 16'^ data the pole fits generally have 
a poor x^/dof. We also find that the simultaneous and 
sequential fits to the and mass dependence for the 16"^ 
data differ at 1.2a. Consequently we only use the 16'^ 
data to check that the finite-volume effects are small. 
Tables Hill and HVl demonstrate that this is the case. 

We present the results from a fit to the 24'^ x 64 data 
sets using ^ in Fig. [1] Here the curve shows the fit func- 
tion at the physical meson masses, while the difference 
/o(g^m^^",TO'^") - /o(g^mP'^y^TO^^y") has been sub- 
tracted from our raw data points and the small scatter 
is indicative of the quality of our fit. 

The quark mass dependence of (O is presented in 
Fig. [2l The solid line represents the fit function eval- 
uated at = 0, plotted as a function of m^, while the 
dashed line is the contribution coming from the 0{p'^) 
terms in the chiral expansion, 1 + f2- Our results clearly 
indicate a sizeable, negative value for A/ = —0.013(3), 
in contrast to the recent ChPT based results of 0, BISl- 
In Fig. [2] we also overlay the results given in Table HVl for 
/o(0) obtained from individual pole fits on each of our 
ensembles and earlier Nf = 2 results [l3| . 

So far, we have assumed a pole dominance behav- 
ior ([nj in our lattice data, whose q^ dependence differs 



marginally at NLO from an expression obtained in ChPT 
7] . In order to estimate the systematic error due to this 
choice, we also present in Table ITVl results for fits to our 
data using a quadratic ansatz 



/o(9^) = 0.0 + aiq^ + a2q^ 



(12) 



together with a chiral extrapolation using (jlOp . A simul- 
taneous fit similar to ([9]) is possible via 



1 + /2 + (m^ 



K-mlf{Ao + Ai+A2{mj^ 



ml)) 



{As + {2Ao + Ai){mj, + mi)) q^ 
Ao + A^iml + ml)) q^ . 



(13) 



The form of this ansatz is motivated by the expression 
obtained in ChPT 7]. We quote the result from a fit to 



TABLE IV: Results for /+(0) using pole dominance (fTTj) and 
quadratic (|12[) fits to each data set, together with the chiral 
extrapolations using (|10p with the 24"^ x 64 data only. The 
final row gives the results for simultaneous and quark mass 
fits (([9]) and (|13p ') using the same data sets. 



Pole 



Quadratic 



am^d 


/+(0) 


xVdof 


/+(0) 


xVdof 


16^ X 32 


0.03 


0.99925(8) 


5.0/3 


0.99938(12) 


4.2/2 


0.02 


0.9951(6) 


13.5/3 


0.9959(9) 


13.0/2 


0.01 


0.9889(26) 


13.9/3 


0.9866(33) 


10.9/2 


24^ X 64 


0.03 


0.9991(2) 


2.1/3 


0.9990(2) 


1.5/2 


0.02 


0.9960(7) 


2.3/3 


0.9962(9) 


1.9/2 


0.01 


0.9841(29) 


10.4/3 


0.9806(39) 


7.7/2 


0.005 


0.9774(35) 


4.0/3 


0.9749(59) 


2.7/2 


chiral 


0.9644(39) 


3.4/2 


0.9622(61) 


5.1/2 


sim. fit 


0.9644(33) 


28.7/16 


0.9610(43) 


26.4/14 
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the 243 X 64 data using ^ in the last row of Table HVl 
where we find that the results of the two fits, ([9]) and 
(Uni), agree within statistical precision and we take the 
difference (0.0034) as an estimate of the systematic error 
in choosing ^ as our preferred ansatz. 

Recently, an alternative parametrization, obtained by 
using analyticity and crossing symmetry, has been pro- 
posed |;21|] . We find that fitting our data using this ansatz 
leads to results that lie within the systematic uncertainty 
of 0.0034 discussed above. 

Our simulations are performed with a strange quark 
mass (ams + amies — 0.043) which is heavier than the 
physical mass {arus + arrives — 0.037). Both ^ and (fT3|) 
are modelled according to ChPT and this mass differ- 
ence is corrected when we insert the physical kaon mass 
to obtain our final result. This correction is accurate in as 
much as our extrapolation model describes our data, and 
any error introduced is included in our estimate of the 
systematic error. Future simulations will include a sec- 
ond valence strange quark mass to decrease the reliance 
on our fit model. 

Finally, since we simulate at a single lattice spacing, we 
are unable to extrapolate to the continuum limit. How- 
ever, leading lattice artefacts with domain wall fermions 
are of ©(o^Aq^^); assuming Aqcd ^ 300 MeV we es- 
timate these to be no larger than « 4% (of 1 — /+). 
A comparison of the pion and kaon decay constants ob- 
tained from our simulations with their physical values 
provides a test for the reliablity of our result. After in- 
cluding the effects to NLO due to chiral logs, we find /tt 
and fx about 4% below experiment , which is consis- 
tent with our estimated scaling error. We will explicitly 
check this for Kg^ decays on our new ensemble which is 
being generated on a finer lattice. Note that our cur- 
rent uncertainty is dominated by statistics and the chiral 
and extrapolations and not by the discretization error. 
Hence our final result is 



1+fo 



/+(0) = 0.9644(33) (34) (14) , 



(14) 



where the first error is statistical, and the second and 
third are estimates of the systematic errors due to our 
choice of parametrization ^ and lattice artefacts, re- 
spectively. To put this result in context, we com- 
pare our value with other determinations of /+(0) in 
Fig. [31 We see that our result agrees very well with 
the Leu twy ler-Roos value [6] and earlier lattice calcu- 
lations 



itwy] 

m 
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IT^ In particular, we note that 

our findings prefer a sizeable, negative value for A/ = 
-0.0129(33) (34) (14), in contrast to recent ChPT based 
phenomenological results [7, 8, 9|. 

Using |Ks/+(0)| = 0.2169(9) from PDG(2006) 



IK.I = 0.2249(9)exp(ll)/+(o) , 



(15) 
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FIG. 3: Comparision with other determinations of /+(0). 



and combined with |14d| = 0.97377(27) Q we find 

\Vudf + Ksf + Kb\^ = I - S, (5 = 0.0012(8), (16) 

compared with the PDG(2006) result, 6 = 0.0008(10). 
Further reduction in the lattice error is imperative. Our 
interpolation systematic is removable in principle (2^ 
and we are in the process of addressing both this and 
discretization systematics with a new set of simulations. 
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4 Using the result from Q gives \Vus\ = 0.2247{5)(11). 
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